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Abstract -The interplay between Coulomb friction and random excitations is studied experi- 
mentally by means of a rotating probe in contact with a stationary granular gas. The granular 
material is independently fluidized by a vertical shaker, acting as a "heat bath" for the Brownian- 
like motion of the probe. Two ball bearings supporting the probe exert nonlinear Coulomb friction 
upon it. The experimental velocity distribution of the probe, autocorrelation function, and power 
spectra are compared with the predictions of a linear Boltzmann equation with friction, which is 
known to simplify in two opposite limits: at high collision frequency, it is mapped to a Fokker- 
Planck equation with nonlinear friction, whereas at low collision frequency, it is described by a 
sequence of independent random kicks followed by friction-induced relaxations. Comparison be- 
tween theory and experiment in these two limits shows good agreement. Deviations are observed 
at very small velocities, where the real bearings are not well modeled by Coulomb friction. 



Introduction. — The role of fluctuations in solid-solid 
interactions with friction has been increasingly studied in 
the last years for different systems and different scales. 
Starting from the 60s, engineers studied the effect of noise 
on dry friction and dry contacts as a way to model the 
stability of buildings under earthquakes; see, e.g., [l}{3]. 
More recently, physicists have started studying similar 
problems, but from a more microscopic point of view, by 
looking at the effects of noise on "small" systems with 
few degrees of freedom, such as those studied, for exam- 
ple, in nanofriction experiment s |4| , particle separation J5J, 
ratchets and granular motors [6)48], as well as droplet dy- 
namics on surfaces [9 - 11 , which involves forces similar to 
dry friction. 

Dry or "Coulomb" friction is a nonlinear dissipative 
force that depends on the relative velocity v between two 
sliding surfaces, which, in the simplest model, takes the 
form Ff ricf o-(v), where cr(v) is the sign of v (and zero when 
v = 0) and -F/Hct is the magnitude of the friction or con- 
tact force. The effect of this force, as is well known, is to 
define a threshold for external forces F ext , such that if the 
surfaces are relatively at rest (v = 0), then the external 
force F ext does not induce motion unless \F ext \ > \Ff rict \. 
This is the basis of stick-slip motion, which can be ren- 
dered all the more complex by the introduction of time- 
dependent forces or external noises. P.-G. de Gennes 12 



showed, in particular, that although a Brownian particle 
affected by dry friction can never come to rest with "pure" 
Gaussian white noise, its diffusive properties are very dif- 
ferent from that of normal Brownian motion with only 
viscous friction. In later studies of de Gennes's and other 
models, it was also shown that some features of stick-slip 
motion remain at the statistical level, e.g., in the way in 
which the correlation function or the power spectrum of 
Brownian motion shows a transition as a function of dry 



friction and external forces; see 13 - 17 



Here we report on the experimental observation of some 
of these properties relating to random motion with dry 
friction. In our experiment, sketched in Fig. [T] a macro- 
scopic rotator has its axis suspended to a couple of bear- 
ings, which are the source of dry friction, and is put in 
random motion by immersing it in a fluidized granular 
gas 18 , which provides the noise source. The motion of 



the rotator is recorded at high frequency, and its statis- 
tical properties are analyzed within a suitable theoretical 
framework, which includes a Boltzmann equation model, 
used in the low collision regime, as well as a Langevin 
model in the high collision regime. These two regimes 
are probed by changing the fluidizing properties of the 
granular gas (e.g., the shaking amplitude) and lead us 
to study two very different phenomenologies associated 
with discontinuous (Poisson-type) noises, on the one hand, 
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Fig. 1: Setup and definition of theoretical quantities. 



and continuous (Gaussian-type) noises, on the other. Our 
results compare well with the theoretical predictions ob- 
tained with the models mentioned above and provide a 
first experimental verification of many of these predictions. 
Overall, they also show that dry friction, which is surely 
relevant in the discontinuous noise regime, has also an im- 
portant role in the continuous noise regime, even though 
the rotator is kept in endless motion. 

Experimental setup. — The granular gas used is 
made of N — 50 spheres of polyoxymethylene (diameter 
d = 6 mm and mass m = 0.15g) in a cylinder of volume 
V « 1.9 X 10 5 mm 3 (and number density n = N/V). The 
cylinder is shaken with a sinusoidal signal at 53 Hz and 
variable amplitude (measured by the maximum rescaled 
acceleration a max /g where g is the gravity acceleration). 
Suspended into the g as, a pawl (also called "rotator") of 
total surface £ — 1.2 x 10 3 mm 2 (height h — 15 mm and 
base perimeter S — Ti/h), mass M = 6.49 g and momen- 
tum of inertia I = 353 g mm 2 rotates around a vertical 
axis attached to two ball bearings. The position of the 
rotator is recorded in time by an angular encoder (Avago 
Technologies). It is convenient to introduce the radius of 
inertia Rj = y/T/M of the rotator. See Fig. [ljfor a sketch 
of the system and the definition of some quantities. 

A close analysis of the dynamics of the rotator shows 
that it is well described by the following equation of mo- 
tion: 

lj = -Acr(w) - 7 a w + r) coU (t) (1) 

where A = Ff r i c t/I = 38pm4s~ 2 is the frictional force 
rescaled by inertia, j a = 6pml s _1 is some viscous damp- 
ing rate related perhaps to air or to other dissipations in 
the bearings, and f] co ll(t) is the random force due to col- 
lisions with the granular gas particles. The granular gas 
itself is stationary and (roughly) homogeneous. 

The velocity distribution of the spheres on the plane 
perpendicular to the rotation axis is obtained by parti- 



procedure) and is fairly approximated by a Gaussian, 

4>{v) ~ e -" 2 /(2« 2 ) 5 ( 2 ) 

where the "thermal" velocity vq has been introduced. 
Small deviations from the Gaussian are observed but are 
neglected for the purpose of this study; see 



20 for details. 



We have changed the maximum acceleration rescaled by 
gravity a maK /g from 6 to 20, finding for v values from 
200 mm 2 s~ 2 to 500 mm 2 s" 2 . 

The pawl is further characterized by its symmetric 
shape factor (g 2 ) SU rf — 1.51, where (-)surf denotes a uni- 
form average over the surface of the object parallel to the 
rotation axis (see 21 for details). The restitution coeffi- 
cient between the spheres and the pawl has been measured 
to be a « 0.83. It is also useful to introduce the "equipar- 
tition" angular velocity ujq = Vqc/Ri where e = y/jj- 
Note that, because of inelastic collisions and frictional dis- 
sipations, the rotator does not satisfy equipartition and ujq 
is only a useful reference value. 

Boltzmann equation. — Since the packing fraction 
of the system does not exceed 3%, the single-particle prob- 
ability density function (pdf) p{u>, t) of the angular veloc- 
ity of the rotator is expected to be fully described, under 
the assumption of diluteness which guarantees molecular 
chaos, by the following equation (8j|21[|22l: 



d t p(u, t) 
J\p,cj>] 

W{lu'\u) = pS 



d u [( Acr(w) + ~faOj)p(u, t)] + J\p, (f>] 

doj' W{uj\uj')p{uj\ t) - P {uj, t)f c {uj) 



(3a) 
(3b) 



ds 
~S~ 



dv^(v)9[(V(s) - v) -n]x (3c) 



|(V(s) - v) • h\S[u' -u- Aw(fl)], 
[V(s) - v] • n g(s)e 2 



Au(s) = (1 + a)- 



Ri 



l + e 2 g( S y 



(3d) 



where we introduce the rates W(u/|w) for the transi- 
tion u) — > u)', the velocity-dependent collision frequency 
f c (uj) = J duj' W(u>'\ui), the pdf </>(v) for the gas particle 
velocities, and the so-called kinematic constraint in the 
form of Heaviside step function 0[(V — v) • h], which en- 
forces the kinematic condition necessary for impact. We 
also use the following symbols: p — nh, V(s) — ljz x r(s) 
is the linear velocity of the rotator at the point of im- 
pact r(s) parametrized by the curvilinear abscissa s along 
the outer perimeter of the rotator, h(s) is the unit vec- 
tor perpendicular to the surface at that point, and finally 
g(s) — r(s) • i(s)/Ri with t(s) = z x n(s), which is the 
unit vector tangent to the surface at the point of impact. 
We refer to Fig. [l] for a visual explanation of these sym- 
bols. The collision rule is given by Eq. (3d I 21 . Note 
that in our setup, at homogeneous fluidization, we mea- 
sure pS ~ 0.31 mm -1 . 



cle tracking via a fast camera (see 19 for details on the 
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Fig. 2: (a) Two examples of signal uj(t) for different values of 
in the experiment, corresponding to choices of the rescaled 
maximum acceleration o max /g = 4.1 and 20.5 respectively; (b) 
rescaled experimental pdfs of the angular velocity for a range 
of rescaled accelerations going from 4.1 to 21.2. All other pa- 
rameters are in the main text. 



which estimates the ratio between the stopping time ta ~ 
wo/A due to dissipation (dominated by dry friction) and 
the collisional time t c ~ (n£t>o) -1 [j] A transition at /? ~ 1 
is expected between a regime called the rare collision limit 
(RCL) at /3 _1 <C 1, with the rotator at rest most of the 
time, and a regime called the frequent collision limit (FCL) 
at /3 _1 3> 1, with the rotator always in motion, contin- 
uously perturbed by collisions. The difference between 
these two regimes is illustrated in Fig. [2^,. 

The pdfs of the angular velocity obtained experimen- 
tally for values of spanning the RCL and FCL are 
reported in Fig. [2]d. There is a great variability when 
/3 _1 goes from small to large values, i.e., when increasing 
the shaking amplitude and, consequently, the collision fre- 
quency. At large values of /3 _1 , the pdfs rescaled by luq 
tends to superimpose, a sign that ojq becomes the leading 



1 Talbot et. al. 8] consider a different parameter, namely, T* = 
AI /(pL 2 m,Vg), and consider a very thin rectangular rotator of length 
L and a two-dimensional projection of the system with density p, so 
that our nS is their p(2L), while our Rj is their L/(2\/3), leading 
to the correspondence f3 — > \Z6neT* . 
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Fig. 3: Pdf of the pawl's angular velocity in the rare collisions 
limit (RCL), obtained with a maximum rescaled acceleration 
of the shaker given by a max /g = 6.5 (/3 _1 = 1.5). All other 
parameters of the experiment are given in the main text. The 
theoretical prediction |7| is displayed as the dashed green line, 
where only a is fitted with the experimental data. For refer- 
ence, we also show the prediction of the theory in the diffusive 
limit with and without dry friction as the dashed red and blue 
curves, respectively. 



velocity scale. In order to make a more detailed contact 
with the theory and understand the basic properties of the 
velocity pdf, we discuss next the RCL and FCL regimes 
separately. 

Rare collision limit. — As seen above, the pawl in 
the RCL <C 1) is often at rest, resulting in a peak 

around u — in the angular velocity pdf. To describe this 
peak, we approximate the expected stationary pdf as 

p(oj) = aS(uj) + (1 - a)p smoot h{u), (5) 

where a is a suitable weight, decreasing as grows, and 
P smoothie) represents the smooth part of the pdf. 

This form of stationary pdf has been studied in [8|[l7l 



23 . In the RCL, the dynamics is reduced to independent 
collisions followed by friction-induced relaxations. More 
precisely, at a collision time t the rotator velocity changes 
from to w*, depending on the projected impact velocity 
V = v • n and the projected impact point g = r ■ t/Rj, 
and then relaxes according to u = — Ao~(uj*) until a time 
t + t such that ui(t + r) = 0. In this case, the stationary 
average of any function y(u), restricted to the times where 
uj(t) 0, can be written as 



(y) = pS 



dS 
S 



dv\v\4>{v) / dty[w{t)\. (6) 
Jo 



With this formula, we can calculate the characteristic 
function of uj by taking y = e %ktJJ and then invert the 
transform to retrieve the smooth part p S mooth{u). For 
the Gaussian <j>(v) of variance Vq and the particular shape 
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of our rotator, this yields 
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This result is compared with our experimental data in 
Fig. [3J It can be seen that the agreement of the tail is 
very good, considering that we only fit the overall scaling 
factor a representing the weight relative to the S(u>) contri- 
bution. On the contrary, the central part is not well repro- 
duced. We suspect that the discrepancy at low velocities 
is due to a failure of the Coulomb friction model at those 
regimes. A close inspection of single trajectories indeed 
reveals that the free relaxing rotator frequently comes to 
rest with spurious oscillations, likely to be due to the ball 
dynamics inside the bearings. This observation points to 
an interesting application of studying ball bearings under 
random excitation: in our case, the macroscopic obser- 
vation of p(oj) magnifies microscopic features around the 
zero velocity which would be hard to characterize and un- 
derstand otherwise. 

From the experiment, we can also evaluate the autocor- 
relation C{t) — (w(i)w(O)) and power spectrum, 



S(f) 



E 

3 



(10) 



These are shown in Fig. [5] as red curves. To our knowl- 
edge, no theory is available for these quantities in the RCL. 
We notice that the large frequency decay of the spectrum 
S(f) ~ f~ 2 is compatible with a small time exponential 
decay, while the part at small / deviates from it, sug- 
gesting a more rapid decay at large times. These features 
are recovered in the graph of C(t). Note that the power 
spectrum also shows one of the higher harmonics of the 
shaker frequency (3 x 53 = 159 Hz): it emerges only when 
the main signal due to the dynamics of the rotator under 
the collisions becomes weak enough and disappears in the 
FCL where the energy injected by the collisions is larger. 

Frequent collisions and large rotator mass: 
equivalence with continuous white noise. In the 

FCL, it is useful to exploit the difference of mass (here 
e = \JmjM = 0.15) by taking a further eC 1 limit. Such 
a limit is often called a "diffusive limit" and allows us to 
expand the Boltzmann equation ([3| to obtain a Fokker- 
Planck equation or, equivalently, a Langevin equation for 
the pawl velocity [21] having the form 
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Fig. 4: Pdf of the pawl's angular velocity in the frequent colli- 
sions regime (FCL) obtained when the maximum rescaled ac- 
celeration of the shaker is a max /<; = 11.74. The prediction of 
the diffusive limit theory, Eq. ( |13[ ), is displayed as the dashed 
red line with no fitting parameters. For reference, we also show 
with the dashed blue line the prediction of the diffusive limit 
without dry friction. 



where 7 = 7 a + 7 g , 7 9 represents a granular viscosity, j a 
an air viscosity, T g a granular velocity diffusion coefficient, 
and £ is a Gaussian white noise with unit variance. In our 



setting, 7 ff and T g are given by 21 



lg = (1 + Oi) 



T g = (1 + a) lg 



(12a) 
(12b) 



The study of the Langevin equation (111 was initiated 
in [l], received strong impulse by de Gennes in |12| , and 
was completed in 13 -16 . Its stationary velocity distribu- 



tion reads 24 



p(uj) = A/"exp 



(M+A/7) 2 

r 9 /7 



Af- 1 = JvTgh erfc(A/ v ^T g ). 



(13) 
(14) 



Note that in the limit j3 — > 0, e.g., when dry friction dis- 
appears (A — > 0), the stationary pdf goes to a Gaussian of 
variance T g /j. Moreover, assuming also 7 a <ti 7 g (which 
is consistent with the FCL), one has r g /7 m Wq(1 + a)/2, 
so that equipartition with the gas (uj 2 ) = ojq is satisfied in 
the ideal elastic case a = 1 25 1. 

In Fig.|4j we find good agreement between the pdf above 
and the experimental data. This comparison, obtained 
with no fitting parameters, represents one of the first 



known experimental verification of the velocity pdf (13), 
as well as one of the few experimental applications of 
the diffusive limit of granular kinetic theory. The cusp 
at uj = predicted by the theory is a striking conse- 
quence of the presence of Coulomb friction, and is well 
reproduced in the experimental data. At large velocities, 
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w3>A/7«6s 1 , the pdf recovers Gaussian tails due to 
the dominance of linear viscosity. 



In the diffusive limit, where Eq. (Ill holds, theoretical 



expectations for the autocorrelation of the angular velocity 
and for the power spectrum have been obtained analyti- 
cally in 13 16 . These expressions involve the eigenvalues 



of the Fokker-Planck operator, and are too long to be re- 
produced here. A verification of these expressions, shown 
in Fig. [5j indicates that the Langevin model of Eq. (11) 
offers a good description of our experiment in the FCL. 

Relevance of friction in the FCL. A more re- 
fined comparison of experiment and theory is obtained in 
Fig. |6j where three main observables are plotted against 
the rescaled peak of the velocity distribution, the 
variance of the distribution, and the correlation time r c 
obtained by fitting C(t) ~ exp(— f/r c ). The three figures 
also show the predictions, for each choice of the param- 
eters, of the diffusive limit theory, as given by Eq. (Ill, 



together with the predictions of the same theory in the 
absence of dry friction (A = 0). This last comparison is 
useful to evaluate the relevance of the dry friction term. 

In Fig. [6^,, we display the peak of the experimental pdf 
(again rescaled by u)q) as a function of This infor- 

mation, in the RCL when the rotator is most of the time 
at rest, is an indirect probe of a, since the experimen- 
tal value p(0) includes also (1 — a)p smoot h{0)- This figure 
clearly shows the decrease of the peak as /3 _1 increases. 
Moreover, it shows that, when f3~ 1 > 1, such a peak gets 
closer to the values analytically computed in the diffusive 
limit where A is negligible and the pdf tends to a Gaus- 
sian. Interestingly, the variance of the distribution (w 2 ), 
whose formula in the diffusive limit reads 



(u 2 ) = 



Ar„ 



2 7 7 



2 eX P 



A2_ 



(15) 



adheres to the prediction of the diffusive theory even in 
the RCL, while at large /3 _1 both the experiment and the 
diffusive theory go toward the linear limit, where A is 
negligible. The small discrepancies are likely to be due to 
the finiteness of the mass ratio e in the experiment. 

A different scenario is observed for the correlation time 
r c , as the experiment shows a non- monotonous behavior 
when plotted against f3~ , with r c growing when moving 
from the RCL to the FCL, up to a regime where the exper- 
iment is well described by the diffusive theory (red line). 
This regime displays a slight decrease of r c as fi^ 1 is fur- 
ther increased. It is interesting to note that, even at the 
highest values of /8 , the prediction of the linear theory 
without friction (A = 0) overestimates by roughly 40% 
the experimental results. The correlation time is there- 
fore very sensitive to the presence of the dry friction and 
signals it even when the "static" information coming from 
the velocity pdf is basically not affected. 

Conclusions. — We have discussed in this letter the 
results of an experiment in which a rotator is submitted 




iff 
10"' 



co 10 r Experiment (fJ =7.2) 



I O Q Experiment =0.2) 
— - Theory (diffusional limit) 



10 



10 



f [1/s] 



Fig. 5: (a) Autocorrelation and (b) power spectrum of the 
pawl's angular velocity in the FCL (black) and RCL (red) ob- 
tained with a maximum rescaled acceleration of the shaker cor- 
responding to a ma x/g = 4.1 and 17, respectively. The diffusive 
limit prediction is superimposed as the dashed green line. 



to dry friction and collisions with a granular gas. By tun- 
ing the shaking amplitude at constant frequency, we have 
explored different random, Brownian-like dynamics which 
are either dominated by friction or by collisions, as well 
as the crossover between these two extreme regimes. In 
the rare collision regime (RCL), our data for the veloc- 
ity pdf display a macroscopic fraction of events at rest 
(lu = 0) and non-Gaussian tails at high velocity, which 
are both well reproduced by a Boltzmann collision model. 
In the frequent collision regime (FCL), our results for the 
velocity pdf, autocorrelation and power spectrum are well 
explained by a Langevin model with dry and viscous fric- 
tions 



13 , which can be derived from the Boltzmann model 



in the diffusive limit. In this limit, dry friction tends to be- 
come negligible compared to the other forces; this explains 
why, at very high collision frequency, one recovers a phe- 
nomenology partly explained by the Ornstein-Uhlenbeck 
model, if we exclude the correlation decay. 

To conclude, we remark that our experiment suggests a 
useful way to estimate parameters which are somewhat dif- 
ficult to measure directly, in analogy with Einstein's the- 
ory of Brownian motion which gives access to Avogadro's 
number through a macroscopic measurement. In our case, 



Eq. (15 1 together with the expressions shown in (12) may 



be used to obtain estimates of a or A knowing the other 
parameters and the value of (w 2 ). At the same time, our 
experiment offers a positive test of the ability of kinetic 
theory to predict macroscopic Brownian coefficients, such 
as the viscosity 7 S and the noise amplitude T g . Ongoing 
extensions of our study include the coupling of the system 
with a motor 26 in order to apply an external force and 



investigate its interplay with friction and collisions 13 



and the possibility of observing ratchet effects for asym- 
metric rotators [8 20 27 . 
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